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Abstract 

We relate </>(x, s), the average number of sites at a transverse distance x in the 
directed animals with s sites in d transverse dimensions, to the two-point correlation 
function of a lattice gas with nearest neighbor exclusion in d dimensions. For large 
s, 4>(x, s) has the scaling form -^-/(|x|/i? s ), where R s is the root mean square radius 
of gyration of animals of s sites. We determine the exact scaling function for d = 1 
to be f(r) = ^|erfc(r/\/3). We also show that 0(x = 0,s) can be determined in 
terms of the animals number generating function of the directed animals. 

The directed animals(DA) problem describes the large-scale geometrical structure of 
subcritical directed percolation clusters, and qualitatively models diverse situations such 
as trees, river networks and dilute polymers in a flowing solvent [1, 2]. Formally, it corre- 
sponds to p — > limit of the directed percolation problem [3] . The imposed directionality 
changes the universality class, and the large-scale structure of directed animals is different 
from isotropic animals. The enumeration of directed site animals in d + 1 dimensions is 
related to hard-core lattice gas(HCLG) at negative activity with repulsive interactions in 
d dimensions and the Yang-Lee edge problem in d dimensions [4, 5, 6, 7]. This equivalence 
helps to obtain several exact results for the directed animals problem [8, 9, 10, 11]. In 
particular, the critical exponents 9 and the transverse size exponent v± are known exactly 
in d = 0, 1 and 2 [8]. For reviews of available exact results, see [12]. The problem is also 
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Figure 1: Schematic figure of a Directed Animal of size 20 on a square lattice 

related to a recently studied model of quantum gravity [13], and there is an unexpected 
relation between the number of distinct eigenvalues for Potts model partition function on 
strips of width w and the number of directed animals with w sites in two dimensions [14]. 

In this paper, we extend the known relation between the directed site animals enumera- 
tion (DSAE) problem to hard-core lattice gases at negative activity to obtain the average 
number of sites at a given transverse distance x from the origin for d + 1-dimensional 
directed animals from the density- density correlation function of the lattice gas in d di- 
mensions. For d—1, using the exact generating function for this correlation function, we 
determine the scaling form for the average number of sites at a given transverse distance 
in a 2-dimensional directed animal having s sites, for large s. In this case, the average 
transverse size of the animal scales as s^ x , where v± is known to be \ for d — 1. For large 
s, the average value of q-th transverse moment < | a; | ^ > varies as R q s C q , where R s is the 
root mean square radius of gyration of animals of s sites. Using the exact scaling function 
we are able to determine the universal constants C q for all q. 

We consider directed animals on a square lattice. Each site (x, t) has two bonds 
directed outwards towards the sites (x + 1, t + 1) and (x — \,t + 1) (Fig. 1). A directed 
site animal A rooted at the origin is a set of occupied sites including origin, such that for 
each occupied site (x,t) other than the origin, at least one of the two sites (x — \,t — 1) 
and (x + 1, t — 1) is also occupied. Fig. 1 shows a directed animal on a square lattice in 
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two dimensions. The number of sites in A will be denoted by |.4.| or s. We define n(x\A) 
as the number of sites of A having the transverse coordinate x. In Fig. 1, n(x\A) takes 
the values 1,3 and for x equal to —4,2,5 respectively, with |^4| = 20. We shall define 
the radius of gyration of A as x 2 , where the sum is over all sites of the animal. The 

animal shown in Fig. 1, has a squared radius of gyration 21/4. 

We define A(y) as y'- 4 ' summed over all animals A rooted at the origin. If A s is the 
number of distinct animals having s sites, it is easy to see that 

oo 

A{y) = Y.y^=Y. A sV s (1) 

A s=l 

For large s, A s varies as \ s s~ e , where 9 is a critical exponent. The radius of gyration 
R s is expected to vary as s v± , where the exponent u± is related to the animals number 
exponent 6 by the hyperscaling relation 6 = dv± [15]. 

For a directed animal of size s we define the generating function \I/(x, y) as 

vJ/(x;y) = ]Tn(x|.4)^l (2) 

A 

where the summation over A is the summation over all animal configurations. This can 
be written as 

tf(x;y)=5>(x, S )A s y s (3) 

s 

where 0(x, s) is the value of n(x, A) averaged over all animals A of size s. 

The DSAE problem in d + 1-dimensions is related to time development of thermal 
relaxation of a hard core lattice gas (HCLG) with nearest neighbor exclusion on d dimen- 
sional lattice [5] . On a d dimensional body-centered hypercubical lattice (a d dimensional 
body-centered hypercubical lattice is a hypercubical lattice having 2 d ~ 1 bonds going out- 
ward from each site), the dynamics of the lattice gas is defined as follows: The evolution 
is stochastic discrete -time Markovian. At odd(even) times, all the odd(even) sites are 
examined in parallel, and if a site has all neighbors empty, its occupation number is set 
to 1 with probability p, and to with probability (1 —p). The rates of transitions satisfy 
detailed balance condition corresponding to the Hamiltonian 
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H = + 00 U i n j ~ ( ln2 H U i ( 4 ) 

<ij> i 

with z — p/(l — p). 

Let p{p) be the average density of particles in the steady state of this system. In [5], 
it was shown that we have 

My) = -p(p = -v) ( 5 ) 

For the d — 1 nearest-neighbor-exclusion lattice gas, it is straight forward to determine 
the average density corresponding to chemical activity p/(l —p). This for square lattice 
gives 

1 



m = - 2 



1 



1 + y 



(6) 



The animal numbers A s , which are the coefficients in the Taylor expansion of A(y) for 
square lattice can be written as [8] 

A s = / ^{l + cosd){l + 2cose) s -\ (7) 
Jo 2n 

For large s, A s ~ -^=3 s s~^, which shows that in this case A = 3, and 9 = 1/2. 

The derivation of Eq.(5) is easily generalized to the case where the values of p at 
different sites are different. Let the probability that site % gets occupied at time t + 1 
given that all its neighbors are unoccupied at time t be pj. Then the rates of this process 
still satisfy detailed balance condition corresponding to the hamiltonian 

H = +00 Y n i n j - ^2Q nz i) n i ( 8 ) 

<ij> i 

where Zi = Pi/(1 —pi)- The probability that site % is occupied in the steady state depends 
on the p/s for all sites j, and will be denoted by p%({pj})- In the corresponding DA 
problem, we have to define the weight of an animal A as product of weights of all occupied 
sites, the weight corresponding to a site with x-coordinate j being yj. Then define Ai({yj}) 
be the sum of weights of all animals rooted at i. Then clearly Ai({yj}) is a formal power 
series in the variables {%}. If all yj = y, this becomes independent of i, and reduces to 
the function A(y). For unequal y/s, Eq.(5) becomes 
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MM) = -PidPj = ~Vj}) (9) 

Applying the operator y x ^- on the weight of any particular animal gives us the weight 
multiplied by the number of occupied sites with transverse coordinate x in the animal. 
Thus, clearly, we get 

v*-^M{yj})\{v i =v} = *(x;v) (10) 

Let Q({zj}) be the grand partition function for the HCLG given by the Hamiltonian 
in Eq. (8) with Zj = Pj/ (1 —pj). Then, Q({zj}) is a linear function of each of the variables 
Zj. Let rjj be the indicator variable taking value 1 if the site j is occupied, and if not. 
The density-density correlation function of the gas G(i, k) is defined as 



d 

G(i,k; {zj}) =< rjkiji > - < i] k >< r]i >= z k —pi({zj}). 



When Zj = z for all j, the correlation function depends only on (i — k) and hence can 
be written as G(i — k; z). 
Now, using Eq.(9), we get 



tf(x;y) 



1 



-G x; z 



y 



(12) 



i + y V l + y, 

This equation is the main result of this paper, and equates the pair correlation function 
(j(x; z) of the HCLG with the function \&(x;?/) which gives the density profile of the DA 
problem. 

For the special case d — 1, it is a simple excercise to calculate G(x; z) explicitly, using 
transfer matrix methods [16]. This gives 

-i 1*1 

II - — • / 1 ! ■■ 

G{x-z) = 



l + Az 



1 - ^l^Tz 



i + v / TT4i 



(13) 



and we get the explicit expression for ty(x; y) on square lattice to be 



y 



(l + y)(l-3y) 









[ V i + 







-\x\ 



(14) 



This determines the density profile in the constant fugacity ensemble, where an animal 
having s sites has weight y s . However, it is more instructive to look at the profile in the 
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constant-s ensemble. This is obtained by looking at the Taylor coefficient of y s in the 
above equation. This can become rather messy. However, the behavior for large s can be 
determined easily. 

In general, for any dimension d , for large s, 0(x, s) has the scaling form 

(15) 

Usually the argument of the scaling function is defined only upto a multiplicative 
constant. We have made specific choice for this by using the variable as |x|/i? s . The 
normalization of scaling function /(r) is chosen such that is satisfies 



/oo 
= 1 (16) 
-oo 

/oo 
d d x |x| 2 /(|x|) = 1 (17) 
-oo 

In the DA problem, as the number of animals grow as A s , the series expansions for 
\l/(x; y) or A(y) in powers of y converge for y < y c = 1/A. For y near y c , the singular part 
of the function A(y) varies as (1 — yX) 8 ^ 1 . For the HCLG problem, this corresponds to 
a singularity in the series for the density p in powers of activity for the activity zly = 
— (1 + A) -1 [Eq. (12)]. This singularity on the negative real line is the Lee- Yang edge 
singularity for this problem. 

For z near the critical value z = zly, this correlation function G(x; z) is expected to 
have the scaling form 

G(x; z = ZLYe~ e ) = c £~ a g{b\x\e u ) + higher order terms in e (18) 

where b and c are non- universal, lattice dependent constants. We choose b such that 
#(0 = ex P(~0 f° r large £, and c is fixed by requiring g(0) = 1. 

The scaling function tends to a constant limiting value as £ tends to zero, and 
decreases to zero exponentailly fast as £ tends to infinity. Let the power-series expansion 
of g(£) about £ = be given by 

oo 

9(0 = T,9 k e (19) 
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Substituting this in Eq.(18), we get 



k=0 



^ \ uk—a 



ZLY 



G{-K-z)~Y.9kb k \*\ k 1 (20) 



To get 0(x, s), we need to determine the coefficient of z s in the above expansion. From 
the binomial expansion of (1 — -^) vk ' a we immediately get 



G(x;l) ~ 5 fe) £^ |x| r(/ + i)r(a-^) (21) 



For fixed k, for large s, we have 



r( s + i) 

Using this in Eq.(21), we get 



T(s + a- vk) ^ s „_ 1 _„ t (22) 



where the correlation length exponent z/ for the HCLG problem is the same as the trans- 
verse size exponent u± for directed animals. 

For d — 1, it is easy to determine z) explicitly. In this case, using Eq.(13) we get 
the scaling form for G(x; z) as 

1 

G(x; e) = — exp(-2xv^) (24) 
so that a — 1, u — 1/2, b — 2 and c = 1/4. The scaling function g(£) is simply given by 

g(0 = exp[-e] (25) 

hence = l/T[k + 1]. Using this, and the values of a and v in Eq.(23) we get the 
leading singular behavior of <f>(x, s)A s on square lattice to be 

-i oo ( \^3% \k -, 

^■«^ = ^£ r (t + i)r(i-*) (26) 

Since T(l — |) has poles when k is an even integer, only the odd terms contribute to 
the sum. The resulting series is easy to sum explicitly, giving 



k 
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1 / \/3x\ 

(f>(x, s)A s = -3 s erfc ( ) correc ^ on t° scaling terms (27) 

where erfc(x) = e~ x ^ dx. 

This gives <p{x, s) for large s to be 

, , . V37TS / \^3x\ . ,. . . 

<p{x, s) = — - — eric — -j= + correction to scaling terms (28) 

From 0(x, s) we can as well derive the expression for gth transverse moment of the 
directed animals. The gth transverse moment of a cluster of size s denoted by /i qs is 
defined as 

/V= E <f>(hs)\i\ q (29) 

i=— s 

Using the scaling form of (f>(i, s), for large s, we see that {ji qs / sR q s ) is a universal constant. 
Denoting it by C q , we get 



qs 



sRi 

For d — 1, we get 



/oo 
^x|x|V(|x|) (30) 
-oo 



1+ § 



(31) 

(/tl L ^ J 

and i? s = f v^- Hence we get the scaling function(16) in d = 1 to be 



/w =^H7i) (32) 

In the entire low-density phase of the HCLG, in any dimension d, we expect the 
correlation function G(x; z) to have an exponential decay at large |x|. But the behavior 
of the scaling function f(r) for large r is in general different. Suppose ln/(r) varies as 
— r a for large r. Putting this behavior in Eq.(3), and using A s ~ X s and Eq. (12) , we get 

G(x; ZLve-') ~ ]T exp(-|x| a s-^ - es) (33) 

For large |x|, the integral can be estimated by steepest descent, and gives logg(£) 
varying as £t+^. Since this should be linear in £, we see that a = j^. As a check, we 
see that in d — 1, v — 1/2, and /(r) varies as exp(— r 2 ) for large r. 
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The case x = is special, in that the density- density correlation function G(0, z) is 
always equal to p(l — p) for hard-core lattice gas for any d dimensional lattice. Hence, if 
one knows p as a function of the activity z ( equivalently, in the DA problem, one knows 
the animal numbers generating function A(y)), then one can determine \l/(0;y) in terms 
of A(y) alone. For ad+l dimensional DA on a body-centered hypercubical lattice it is 
given by 

*(0-y) = -±-A(y)[l + A(y)] (34) 
1 + y 

For d — 1, using the A(y) from Eq. (6), we get for the square lattice DA problem 



Expanding in powers of y we get 



V 4- (— 1 V -1 

0(0, s) = + \ A } (36) 



where A s is as given in Eq.(7). For large s, 0(0, s) varies as s 1//2 as expected. 

Similar analysis can be extended to higher dimensions. In general the scaling function 
f(r) tends to a finite value as r tends to 0, and hence 0(0, s) varies as s 1 ^ 9 . 

For d > 7, mean-field theory becomes asymptotically exact [19] with R s ~ s 1 / 4 and 
the scaling form of /(r) is 

where Qd is the surface area of a d dimensional unit sphere and is equal to Qd = r(d/2) ■ 
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